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Abstract. Following Sun and Moll [4], we study v p (T(N)), the p-adic valuation of the 
counting function of the alternating sign matrices. We find an exact analytic expression for 
it that exhibits the fluctuating behaviour, by means of Fourier coefficients. The method is 
the Mellin-Perron technique, which is familiar in the analysis of the sum-of-digits function 
and related quantities. 



which is famous because of the enumeration of the Alternating Sign Matrices. Their paper 
provides also some historic remarks about the fascinating story of this subject. 
Sun and Moll write: 

Given an interesting sequence of integers, it is a natural question to explore 
the structure of their factorization into primes. This is measured by the 
p-adic valuation of the elements of the sequence. 
Indeed, in [I] the p-adic valuation v p (T(N)) is studied, i.e., v p {m) denotes the maximum 
k such that p k divides m. 

Most results are, however, for p = 2 only. 
Their key result is (Corollary 2.2): 
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1. Introduction 



Sun and Moll [4J consider the counting function 
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where S p (k) is the sum-of-digits function of k to the base p. 

We use an analytic approach (the Mellin- Perron technique and its extension, [21 E]) to 
derive an exact analytic expression for this function, and this works for all primes. The 
periodicities are made fully explicit by computing the relevant Fourier coefficients. 



2. Results 

Now we study the function v p (T(N)) more closely. The prime p is fixed throughout this 
paper, everything may depend on it, including implicit constants in the O-notation. 
We will prove that asymptotically, 

v p (T(N)) = N\og p + iV$(log p iV) + O(VN), (1) 



where 



fcez\{o} 

is a 1-periodic and continuous function of mean with Fourier coefficients 

(l-2 1+ »+3^)C(Xfc) 2A;ttz 

c k = 7— vj with Xk = ; 

Xfc(l + Xfc)logp logp 

for k e Z\{0}. 

The O-term in Q can also be computed explicitly. Here, the cases p = 1 (mod 3), 
p = — 1 (mod 3) and p = 3 differ in several aspects, thus we formulate one theorem for 
each case: Theorems [TJ [2] and |3j respectively. 

The occurring quantities are expressed in terms of the Hurwitz ( function defined by 



C(s,a) = 



n + a) s 

n>—a x ' 

where we allow arbitrary a£l. For < a < 1, this coincides with the usual definition. 
In our version, however, it is 1-periodic in the second variable, which will be useful for our 
purposes. The special case a G Z corresponds to the Riemann ( function C( s )- 

Theorem 1. Assume that p = 1 (mod 3). Then the asymptotic expansion (which is also 
exact) of the function v p (T(N)) is given by 

v p {T{N)) = AHog p + iV$(log p iV) + *(N) + ± \og p N + f (N). 

Here, 

• the function 

= c k e 2knix 
fcez\{o} 

is 1-periodic and continuous; 
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the function 

V(N) = [N+ I)Vi(log p (iV+ D) + (N- ^-i(log p (iV- D) 

-A^ (log p AO, 
is expressed in terms of the continuous 1 -periodic functions 

i/>j(x)= d ^ 2kmx , je {-1,0,1} 

kez\{o} 

and satisfies 

q?(N) = 0(N 1/2 log N); 
the least order term is given by 



3 * p \r(l)J 6 91ogp 6(p-l) 

with 

ftW - (l + ^)iVlog p (l + ^) „^ ^ + O(l), (2) 
t/ras fo(N) converges to 



3 a '\r(l)J 6(p-i) 

/or iV — > oo . 

TTie Fourier coefficients of$(x), ip±i(x) and ipo are given by 

(1 -2 1+ » +3»)C(Xfc) 2fc7ri 

Ck = 7— VI ™ = 1 ' 

Xfc(l + Xfc)logp 
, _ , , , _ (3 Xfc -l)C(x,) 

«fc,0 — + «fe -1 



Xfc(l + Xfc)logp' 



respectively. 



The case p = — 1 (mod 2) is similar, but the structure of ^(rr) is more complicated: it is 
now composed of 2-periodic functions instead of 1-periodic functions. On the other hand, 
the least order term is less complicated. 

Theorem 2. Assume that p = —1 (mod 3). Then the asymptotic expansion (which is also 
exact) of the function v p (T(N)) is given by 

v p (T(N)) = N\og p + iV$(log p iV) + *(N) + g^j- 
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• the function 

$(x)= Yl c ^ 2kmx 

fcez\{o} 

is l-periodic and continuous; 

• the function 

*(N) = (n+ ^i(io gp (iv+ i)) + (^v - D^-i (iog^(i>r - 1) 

-A^ (log p A0, 
zs expressed in terms of the continuous 2-periodic functions 

fy(x)= Yl d k,e kmx , j G{-1,0,1} 
fcez\{o} 

and satisfies 

\&(N) = 0(N 1/2 log N). 
The Fourier coefficients of §(x), ip±i{x) and ipo are given by 

(l-2 1+ **+3**)C(x fc ) ... 2km 

Ck = 7— r~j with Xk = ; , 

Xk{l + Xk) logp logp 

. axk/2,jm + (-i) k axk/2,-jm 

4,j = o 7— VT"^ /or j G {±1}, 

, , , / (3Xfc ri; 1K( n /2) , if k is even, 

[0, a/ is odd, 

respectively. 

Finally, for p = 3, the result is that the O-term in ([I]) can be omitted: 
Theorem 3. Then Vs(T(N)) is given by 

v 3 {T(N)) = iv(log 3 2 - i) + iV$(log 3 AT), 



where 



$(x)= ^ c fc e 2fc ™ 



fcGZ\{0} 

is a l-periodic continuous function with Fourier coefficients 

2(l-2^)C(Xfe) ... 2A;ttz 

Xfc(l + Xfc)logp logp 

As examples, Figures [IJ [2] and [3] show v p (T(N))/N (gray dots) and an approximation of 
log p ^ + $(log p N) (black line) for p = 2, p = 3 and p = 7, respectively. For the approx- 
imation, 400 Fourier coefficients of $ are used. Note that for p = 3, the approximation 
is much better as the only approximation consisted in truncating the Fourier series, in 
contrast to the lower order terms present in the other cases. 
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Figure 1. Periodic fluctuation for p = 2 



0.5 



0.4 



0.3 



0.2 



0.1 




3 4 5 6 

Figure 2. Periodic fluctuation for p = 3 



Our analysis builds on a result of Delange pQ: 
Lemma 2.1 (Delange PQ). 



JV-l 



5 » = 1o Sp N + ^ + Ar$(1) dog p iV), 



n=0 



with 



Cn = — (l0g27T - 1) - 



2 log p 4 
a l-periodic continuous function of mean zero 

$a)( X ) = ^4V fc ™, 
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Figure 3. Periodic fluctuation for p = 7 



and 



p — 1 



/or fc 7^ 0. 

This already leads to the simplification 

P 



Xk 



2kiri 
\ogp 



v P (T(N)) = ijVlog p iV + iv(log p 2 + ^ I ) 

i T 2iV$( 1 )(log p 2iV) - J_iV$W(log N) + J—^S p (3j + 1). (3) 
— 1 ^ p — 1 ^ I — v L — ' 



N-l 



P — I p 

In the rest of this paper we will analyse the sum 



1 N ~* 
— 5>(3j + l). 



For p = 3, we have S3(3j + 1) = 1 + Ss(j), so that another application of Lemma 2.1 
directly yields the result described in Theorem |3j 

In the following, we will therefore assume that p ^ 3. However, we will not follow 
Delange's original (elementary) approach, but rather use the Mellin-Perron approach as 
described in [2] and extended in [3J. 



3. Rewriting the Remaining Sum as an Integral 
The aim of this section is to prove the following lemma. 
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Lemma 3.1. We have 



r^X><-->-r^(f-^) + ' 



p ^ Fy ' 1 - p V 2 
with 

" 2+io ° ■ jy+i ds 



s(s + l) 



1 pz+ioo ■ 

and 

Ai(s) = 3 p2 a _i — . (5) 
where u G {±1} is chosen such that p = u (mod 3). 

Proof. We intend to use the Mellin-Perron summation formula in the version 



i i-C+ioo I 
- ,/C-ioo y - 



" (/ " V + < +1 * , (6) 



(n + a) s / v y s(s + l) 

-a<n<N x n>— « 

where a G M and C is in the half-plane of absolute convergence of the Dirichlet series 
J2 n >-a (n+a) s > c ^ K (4-3)]. Note that we do not impose the frequently used restriction 
< a < 1, but we sum over n > —a, which amounts to the same. This version, however, 
smoothes the following calculations. 
For j G {—1,0, 1}, we set 

S p (3n + j) - S p (3n+j - 1) - 1 

% H = ; 

1 — p 

and 

1 „>^ /3 (»+i) 

Note that Aj(s) converges absolutely for Res > 1. The additional summand —1 in the 
numerator of %(n) leads to simpler expressions for Aj-(s), the denominator also simplifies 
the expressions and is also present in 

To prove Q, we apply (|6| with a(n) = aj(n) and a = j/3 for j G { — 1,0, 1}, sum up 
the result and use Abel summation. This yields 

N-l 



I = - n)(ai(n) + a {n) + a_i(n)) + iVai(O) 

n=l 

1 7V_1 1 

= — J2( N - n )( S p(3n + 1) - S p (3n - 2) - 3) + —N(S P (1) - 5 P (0) - 1) 

n=l ^ 
1 / iV " 2 

= — — ( S P (3iV - 2) + S p (3n + 1)(JV - n - (JV - (n + 1))) 
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N 3N 2 ' 



-(N-1)S P (1) + NS P (1) + 



N-2 



i / iV ~ z N 3N 2 \ 

— (S p (3N-2) + Y,S P (3n + l) + S p (l) + --—) 

P ^ n=l ~ ' 



N-l 



= —<£S,{*n + l) + — (--—), 

y n=0 y 

as requested. 

In order to compute J, we need explicit expressions for the Ay It is well-known (and 
easy to see) that for a positive integer m, we have 

S p (m) - S p (m - 1) = 1 - (p - l)up(m), 

which immediately results in 

Oj-(n) = w p (3n + j). 

Thus we get 

v p (3n+j)_^ x - 1 \- \- 1 



a,w= E ^^ = E E irhv-Y. E 



»>^i/3 ( n+ 3) k>l n>-j/3 ( U+ l) fc>l n>-i/3 ( n + 3) 

fp(3n+j)>fc p fe |3n+j 

By definition of it, we have = (mod 3). Thus 3n +j is divisible by p k if and only 
if there is an integer m such that 3n + j = p k (3m + ju k ). This results in 

1 



a 3 w = E E 



fc>l m >_j U fc/3 P V™ + ~3~/ 

We now split the sum over k according to k's parity and get 

Aj(s) = § m S/3 ^F^F + § m £ /3 p (2fc - 1)s (-+f) s 
p 2s - 1 

as claimed in □ 

4. Computing the Asymptotic Main Terms 

The further strategy is now to shift the line of integration to the left. The residues at the 
poles yield the main terms in the asymptotic expansion. We prove the following lemma. 
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Lemma 4.1. We have 

+ (log, N) + ^(N) + A (2) log p N + /f (8) 
1 r-z+ioo ■ +1 rfs 



s(s + l)' 



1 /'-4+ 400 / 

w/iere 

• ty(N) has been defined in Theorems^ and^ respectively and has the properties 
given in these theorems; 

• the quantities f[ 2 ^ and f^ depend on p modulo 3: 

f (2)_ f(2)__JL 1 fO) 

/x _U ' /o " 18 9(1 -p) W 
if p = —1 (mod 3) and 



f(2) _ 1 

"9' 



r r (l)A , 9i(N)-g-i(N) 1 1 1 



(10) 



3 to nr(§)/ 6 91ogp 18 9(1 -p) 

if p = 1 (mod 3) £/ie convergent function gj{N) defined in (J2J); 
• £/ie 1-periodic function 

$( 2 )( X )= ^ c (2) e 2fc^ (11) 

fcez\{o} 

Cfc = ~^T^ vi fork^O 

zs continuous. 

Proof. We shift the line of integration to Res = —1/4, 

/= ss5:y 1|ao A iW(^+i) ^riy+E^ 

^J-j-bo V I s p ie S 

where r s denotes the residue of the integrand in s. Here, the usual growth estimates for 
the ( function on vertical lines (cf. (5j § 13.51]) 



((s,a) = 0(\t\ T{a) hg\t\) with r(<r) = < 



1 < er, 

I- a 1/2 < «7 < 1, 

1/2 0<<7<l/2, 



1/2 -a o < 



10 CLEMENS HEUBERGER AND HELMUT PRODINGER 

have been used, where a = Re s and t = 1ms, as usual. 
The residue at s — 1 equals 

3iV 2 1 



n 



2(1 -p) 9(1 -pY 
The residue at s = is calculated as 

= ~\ N N + N (-l **> fa + 21^ + i) + /l<2 ' log * " + /f) + ^ 

with the quantities f± and defined in ^ and (10). 



Finally, there is a simple pole of the integrand at s = Xk/2 = kiri/\ogp for k £ ^\{0} 
with residue 

C(x fc /2 ! j73) + (-i) fc C(x fc /2,^73) / jy+x fc/2 

rXfc/2 ' ^ 2 Xfc/2 (l + Xfc/2 )logp I + 3/ 

Approximating (N + j/3) 1+Xfc/2 by A^ 1+Xfe / 2 and correcting the error yields 

= Tj=-i (C(x fc /2, j/3) + (-i) fc C(x fc/2 ,^73)) ^ 1+ ^ /9 

TXk/2 ' 2 Xfc /2(l+Xfc/2)l0gp 



C(xfc/ 2i j73) + (-i) fc C(xfc/2,^73) f( N + i\ 1+Xk/2 _ N i+ Xk 

^ 2 Xfc/2 (l + Xfc/2 )logp U 3/ 

As ((s, —1/3) + £(s, 0) + 1/3) = 3 s £(s), the first summand equals 

1 + 3Xfc/2 C(x fc /2) Arl+Xt/ , 

2 Xfc/2(1 + Xfe/2)logp 
which leads to no contribution for odd fc. Summing over all ^ yields 

^r Xfc/2 = iV$( 2 )(log p iV) + vI/(iV), 



/2 



(13) 



where the 1-periodic function <3>( 2 ) is defined in (11) — it corresponds to the first summand 
in (13) — and ^(N) is the function defined in Theorems [l] and [5J 



As 

(N + |) +Xk/2 - N 1+x ^ = 0(mm{N, |Jfe|}) 



and 



((Xk/2,j/3) + (-l) k aXk/2,uj/3) 

2Xfc/ 2 (l + XA:/2)logP 



0(|A;|- 3/2 1oe 



by (12), we obtain 

tf(iV)< \k\~ 1/2 ^g\k\ + N\k\- 3/2 log\k\ 

l<|fe|<JV \k\>N 
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< N 1/2 log N + N 1/2 log N < N 1/2 log N, 

as claimed in Theorems G] and [21 

As all involved Fourier series converge absolutely and uniformly, the periodic functions 
are continuous. 

□ 

5. Remainder Integral 



Finally, we deal with the remainder integral in ([8]). It is clear that R = 0(N 3 ^), but 
even better, it can be computed explicitly 



Lemma 5.1. We have 



1 ^ r-k+i°° / jy+i ds 2p 



,-=_!•' -4-100 



s(s + l) 9(p-l)" 
Proof. Expanding the denominator p 2s — 1 in a geometric series yields 



r- 

■2s 



--iEi/;:(cH) + «(4)>i^D : 



s(s + l) 



j=-l £>0 

where exchanging summation and integration was legitimate due to absolute convergence. 

Shifting the line of integration back to Res = 2 and taking the residues in the simple 
poles at s = and s = 1 into account yields 



J 

We have 

f2+ioc 



3J s(s + r 



1 r ds, v {)p-(N + ^ +i ds 

J2-ioo V 6 7 V 



2?r« A-ioc v ' 3 / V 3/ s(s + l) 

p -m r2+ioo . (p m - V )j , d S 



r c(4)(^ + 



27 ™ A-ioc V ' 3 / V 3 3 / s(s + 1) 

-i)j/3<n<p m Af+(p m -i))j/3 

by the Mellin- Perron summation formula again. In the cases of interest, we have either 
v — u and m odd or v — 1 and m even, which implies that p m — v is divisible by 3 in our 
cases. Thus we obtain 



I C(s,—)p ms (N + , 
2mJ 2 _ ioD H 3 J 1 \ 3J s(s + l 
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p . m {p m N + (p m - v)j/3 + \vj/3] - l)(p m N + (p m - v)j/3 + \vj/3] ) 



2 

Summing over j yields 



r ^ ' 9(p-i; 



£>0 

□ 

Combining O), Lemma 3.1, Lemma 4.1 and 5.1 yields the results given in Theorems [l] 
and 13 
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